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CONVERGENCE IN Lp FOR FEYNMAN PATH INTEGRALS 


FABIO NICOLA 


Abstract. We consider a class of Schrddinger equations with time-dependent 
smooth magnetic and electric potentials having a growth at infinity at most 
linear and quadratic, respectively. We study the convergence in LP with loss of 
derivatives, 1 < p < oo, of the time slicing approximations of the corresponding 
Feynman path integral. The results are completely sharp and hold for long 
time, where no smoothing effect is available. The techniques are based on the 
decomposition and reconstruction of functions and operators with respect to 
certain wave packets in phase space. 


1. Introduction 

Feynman path integrals were introduced in 1948 [T8| fT^ to provide a new formu¬ 
lation of Quantum Mechanics and nowadays represent a fundamental tool in most 
branches of modern Physics. In particular, R. Feynman suggested the construction 
of the integral kernel K{t, s,x,y) of the Schrddinger propagator as a suggestive 
sum-over-histories, in the following sense. First of all the kernel K{t,s,x,y) itself 
is interpreted as the probability amplitude for a particle to be at the point x at 
time t provided it was at y at time s {x,y ^ M'^). Now, in the computation of 
this quantity every path 7 joining y and x, therefore satisfying 7(5) = y, 'j(t) = x, 
carries a contribution which is proportional to where S'[7] is the action 

along the path 7: 

>^[7] = ^ ^( 7 ('r), 7 (r),r)dr, 

L being the Lagrangian of the corresponding classical system. The total amplitude 
is hnally obtained by superposition and can be written symbolically as an integral 

K{t,s,x,y) = j 

over the space of paths satisfying the above boundary conditions, where is 
a suitable measure on this space. Although such a measure actually does not 
exist in the measure theoretic sense [6], several rigorous justihcations have been 
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proposed by many authors and from different viewpoints (analytic continuation 
of the parabolic propagator, inhnite dimensional oscillatory integrals, stochastic 
integrals, etc.). The literature is enormous and we refer to the books [Tl HSlHGlHS] 
and the references therein. Instead here we focus on the original approach of 
Feynman [TSlIT^ via time slicing approximations, which was carried on in a rigorous 
way in the papers [201 - I2T11261 - 12^13^ 135113^411151115^ . Briefly one argues as follows. 
Suppose that for \t — s\ small enough there is only one classical path 7 (i.e. a 
path satisfying the Euler-Lagrange equations) satisfying the boundary condition 
7 ( 5 ) = ?/, 7 (t) = X. Dehne then the action 

(1) S{t,s,x,y)= f L( 7 (r), 7 (r),r)dr, 

J S 

along that path. 

Consider the operator {t, s) dehned by 

( 2 ) s)f(x) = ^ dy. 

{2m{t - s)hp^ J^d 

The idea is that this operator should represent a good approximation of the actual 
propagator when \t — s| is small (in fact, for the free particle E®(t,s) coincides 
with the exact propagator). In general one then considers a subdivision : s = 
to < ti < ... < = t of the interval [s, t] and the composition 

(3) t, s) = tL.,)E^^\tL-i, tL- 2 )... s), 

which has integral kernel 


( 4 ) 


K^^\Q,t,s,x,y) 
n (27ri{tj - 


- 1 ) 


L 

exp (^h~^ S{tj, fj_i, Xj,Xj_i)^ 
i=i 


L-l 

Y[dxj, 

1=1 


with X = xl and y = Xq] see Figure 1. 



Figure 1. A piecewise classical path in spacetime. 
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Feynman stated that the operator s) converges (in some sense) to the 

actnal propagator as a;(r2) := snpj^;^ — tj_i| —)■ 0. In view of the semi- 

classical approximation one also introdnces higher order parametrices 
N = 1,2 .. and the corresponding time slicing approximations s) (see 

Section 2 below), and set the same convergence problem. In the present paper we 
investigate this issue for the uniform topology of bounded operators on L^-based 
Sobolev spaces. 

We consider the Schrodinger equation 


(5) 


ihdtu 


1 

2 


h'^Au -|- V{t, x)u 


where 0 < h < 1 and the potential V{t,x), t ^ M., x E satisfies the following 
hypothesis. 


Assumption (A) V{t,x) is a real function of {t,x) e R x with dfV{t,x) 
continuous in {t, x) G R x R'^, for every a G and satisfying 

\dfV{t,x)\ < Ca, |a| > 2, {t,x) G R X R'^. 


Under this hypothesis the exact propagator U{t,s) was constructed in [2011^ : 
it turns out that, for \t — s| small enough, U{t, s) is an oscillatory integral operator 
of the form 

(6) U{t,s)f{x) = ^ [ E^~'^^^’'''^'y'>b{h,t,s,x,y)f{y)dy 

[zni[t — s)n) ' j^d 

for some amplitude b such that dfd^b{h, t, s, x, y) is bounded, for every a, (3 E 
0 < h < 1. Instead, for \t — s\ large, U{t, s) can be written as the composition of a 
hnite number of oscillatory integral operators as above. 

Moreover t, s) (for fixed N) was also shown to converge in a quite strong 

topology to the actual propagator as a;(r 2 ) —)■ 0 when t — s is small enough, which 
implies the convergence in the uniform topology of bounded operators on L^(R'^). 
While we are not aware of similar results in for p 7 ^ 2, it is clear that in that 

case the propagator is in general not even bounded on but a loss of derivatives 
occurs. For example, for the free propagator we have 

||e*^*^/||LP < C'lKl - hA)’^/‘^f\\LP, k = 2d\l/2 - l/p\, 1 < p < 00 . 

This estimate was proved in |31] in the case h = 1, whereas the general case follows 
at once by a scaling argument. Generalizations to different classes of potentials 
were proved, except for the endpoint, in [ll[7l[3ni[3T],|l5] , exploiting in an essential 
way some smoothing effect (Gaussian estimates for the heat kernel). The above 
loss of derivatives is optimal [5] and is explained by the fact that the characteris¬ 
tic manifold of the Schrodinger operator, that is a paraboloid, has non-vanishing 
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Gaussian curvature. For hyperbolic equations [I71ll7lll9] we have instead the loss 
k = 2(d — l)|l/2 — l/p|, because there is one flat direction in that case. 

Motivated by this model, we will focus on the convergence of the time slicing 
approximations in the following Sobolev spaces. 

For 1 <p<cx),A;GM, define the space 

(7) il = = {/ £ S'(R'') : ||/||£j = ||(1 - < «d}. 

Hence, is the usual Sobolev space as a vector space, but with a rescaled norm 
at the Planck scale. Here is our main result. 

Theorem 1.1. Assume the condition in Assumption (A). Let 1 < p < oo, k = 
kp = 2d\l/2-l/p\. 

a) For every T > 0 there exists a constant C{T) > 0 such that, for all f G 

|s — t| < T, 0 < h < 1, 

(8) \\U{t,s)fh.<C{T)\\fh., l<p<2, 

(9) \\U{t,s)f\k._^<C{T)\\fh., 2<p<oc. 

b) For every T > 0, N = 0,1,2,..there exists a constant C{T) > 0 such that, 
for 0 < t — s < T and any subdivision G of the interval [s,f], / G 

0 < h <1, we have 

(10) - C7(i,s))/lli- < C{T)h’^u{ap*'\t - s|||/|(j,, 1 < p < 2, 

(11) ||(B<"l(fi,«,a)-C/(i,s))/||p,^ <C(r)/iVsl)«'"‘l<-s|ll/IU£ 2<p<co. 

We will show in the last section that these are, in general, all possible L^- 
estimates for Schrodinger propagators of the above type. Extensions to the case of 
magnetic fields (cf. [211EIIE2]) or even rough potentials will be discussed below. 

Let us now say a few words about the strategy of the proof. There are two main 
issues. First, for long time the propagator does not have the integral representation 
(]6|) anymore and in general no smoothing effect is available. For example, for the 
quadratic potential V{t,x) = \x\‘^/2 (and h = 1 for simplicity), the propagator 
U{t,0) at time t = kn, /c G Z, has distribution kernel — y) (see 

e.g. [33]). Secondly, the space of bounded operators Lf —?• is of course not 
an algebra, and hence it does not behave nicely with respect to compositions as 
those in fl20D . To overcome these problems we lift the whole analysis to the phase 
space, using ideas originated by the seminal papers [T6l[38lll2],[50] and developed 
by the author et al. in [8l[Tni[T2|. The key tool is a space of operators highly sparse 
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with respect to certain wave packet decompositions. To be precise, consider a fixed 
Schwartz function g ^ 0, and the wave packet 

T^{z)g{y) = e^^^giy -x), z= x 

which is highly concentrated, in phase space, near (a:, We will consider operators 
T satisfying (when h= 1) 

\{TTi{z)g,T^{w)g)\ < (7^(1 + \w - x{z)\)~'^, z,w e 

for every m G N and for some canonical transformation y : —)■ Com¬ 

position of operators corresponds to the composition of the corresponding maps 
x’s. Moreover these operators are trivially bounded on the so-called modulation 
spaces MP, 1 < p < oo (see [igES] and Section 3 below), which measure the phase 
space concentration of a function. We finally come back to spaces by means of 
a nontrivial embedding recently obtained in |37] . 

While endpoint continuity results for p ^ 2 typically require a painful analysis 
in some version of the Hardy space (see e.g. [I61I111I171119]), one of the novelties 
of this approach is that such technicalities are hidden under the above mentioned 
embedding and the analysis is focused on how the operators move the wave packets 
in phase space, which is much smoother and closer to intuition. For this reason 
and in view of the results in [HllISlIlS], we believe that the present approach could 
also extend to the case of rough potentials (and degenerate operators). This would 
provide an interesting low regularity version of the results of |211I52]- We plan to 
investigate these issues in a subsequent paper (see also the last section for some 
partial result in this connection). 

Briefly, the paper is organized as follows. In Section 2 we recall the short-time 
analysis of [2011211152]. Section 3 is devoted to the phase space methods used in 
the sequel and, in particular, to the class of operators, closed under composition, 
mentioned above. In Section 4 we will show that the short-time propagators U{t, s) 
and approximations s) belong to such operator class. Theorem 11.11 will 

be proved in Section 5. Section 6 treats a generalization concerning the convergence 
of s) with iV > 1 in the presence of a magnetic field satisfying the same 

assumptions as in [52] • Finally in Section 7 we show the sharpness of the results and 
provide an extension of the first part of Theorem II.II to a class of rough potentials. 

2. The Schrodinger flow: short time analysis 

In this section we recall the main results obtained in [201 [2T] : see also [52]. We 
will always assume the hypothesis in Assumption (A) in Introduction. 

First of all we fix the following notation. We denote by S'qq the class of smooth 
function a(x, y) in which are bounded together with their derivatives, endowed 
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with the seminorms 

||a||m = sup ||5“5^a||£,oo(R2d), meN. 

\a\ + \l3\<m 

Consider now the Hamiltonian 

+ V{t,x). 

Denote by {x(t, s,y,ri),^{t, s,y,ri)) {s,t G M., y,ri E M'^), the solution of the corre¬ 
sponding system 

X = ^, i = -VxV{t,x) 

with initial condition at time t = s given by x{s, s, y, rf) = y, ^(s, s, y, rj) = rj. The 
flow 


(12) {x{t, s, y, s, y, y)) = x(f, s){y, y) 

dehnes a smooth canonical transformation ^(t, s) : satisfying for every 

T > 0 the estimates 


(13) \d^d^x{t,s,y,y)\ + \d^d^^{t,s,y,y)\<Ca,yiT) y,y eW^, 

for some constant C(T) >0, if|f — s| <T (see Proposition 1.1]). 

Moreover, there exists h > 0 such that for |t — s| <5 and every x,y E 
there exists only one solution such that x{s) = y, x{t) = x. By computing the 
action along this path, as in ([1]), we dehne the generating function S{t,s,x,y) for 
|f — s| <6. It satishes the estimates 


(14) 
and 

(15) 


\t-s\ \d^d^S{t,s,x,y)\ < 


l‘ - *1 -j 


dy^ 


Q^l T 1/^1 ^ 2 


for some 5 > 0 and every x,y E always for |f — s| < 6. 

We now come to the construction of the parametrices. For |t — s| < 6, dehne the 
operator s) for |f — s| < 5 as in ([2]). 

For = 1, 2,... we dehne s) as 

(16) E^^\t,s)f{x) = ^ [ E'^~'^^^’''’^’y^eN{h,t,s,x,y)f{y)dy, 

[zni[t — s)n) ' j^d 

with eN{h,t,s,x,y) = J2k=ii'^^)^~^^k{t,s,x,y), where the amplitudes satisfy 
the transport equations 


dcLk OS OcLk 1 / A ri 

J=l J J 




2 ^x^k—1') 
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with ao(t, s, x,y) = 0 and the initial conditions ai(s, s, x,y) = 1 and ak{s, s, x, y) = 
0 for A; = 1, 2,... 

We have, for every m G M, 

(17) \\ak{t,s,-,-)\\m<Cm for|f-s|<h. 

The operators (t, s) are parametrices in the sense that, for = 0,1,..., 

{zhd, + - V{t,x))E(^\t,s)f = G^^\t,s)f 


with 

( 18 ) GW(t.»)/- ^ f a:,!/)/(!/)da, 

[zm[t — s)n) ' jRd 

where satisfies the estimates Propositions 1.5, 1.6] 

(19) \\gN{h,t,S, •, Ollm < Cn,h^^^\t - 

Moreover, for a snbdivision Q : s = to < ti < ... < ti = t we define 

(20) EW (fi, t, s) = E^^^ {t, tL_,)E^^^{tL-l , tL- 2 ) . . . E^^'^{h , s ). 

As already observed in the introduction for |t — s| < 6 the propagator has the form 

pi) L 

for an amplitude h such that d^d^b{h, t, s, x, y) is of class in t, s, x, y and satis¬ 
fying 

(22) \\b{h,t,s,-,-)\\m<Cm 

for |f — s| < h, 0 < h < 1, m G N. 

When \t — s| is large, U{t, s) is the composition of a hnite number of such oscil¬ 
latory integral operators. 


3. Phase space methods 

We recall here some basic facts of phase space analysis used in the sequel; see [25] 
and especially [I5] for a Mathematical Physics perspective. Then we report on some 
results from p fTOllSO] . 
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3.1. Phase space concepts. For a point z = (x, G x and a function / 
on we denote the phase-space shifts by 

d 

'^{z)f{y) = - x), where = y ■ ^ = ^ y^^^ . 

k=l 

The short-time Fourier transform (STFT), or Bargmann transform, or FBI 
transform, of a function/distribution / on with respect to a Schwartz window 
function g G \ {0} is defined by 


(23) 


L/(2) = (/.>r(s)9) = f f{v)g{v-x)e dy, 

jR'i 


for z = {x,^) X Its adjoint is given by 

V:F{y)= [ F{z)7r{z)g{y)dz, yE 


and we have the inversion formula (see e.g. [25l Corollary 3.2.3]) 

V = —1 


(24) 


/ = 


(2t) 


L2 


(27r)^ 


2 

L2 


y9f{z)T^{z)gdz. 


For 1 < p < cx) the modulation space is the space of distributions 

/ G such that their STFTs belong to the space with norm 


II/IImp(R‘*) ll^9/llLP(]R2d). 

This definition does not depend on the choice of the window g G S{W^),g ^ 0, 
and different windows yield equivalent norms on [251 Theorem 11.3.7]. We also 
have 

Vg : ^ V* : ^ MP(R'^) 

as bounded operators. 

The following embedding results (see m Theorems 1.3, 1.4] and also [36]) are 
crucial in the sequel. 

For 1 < p < cx), fc G R, let = T^(R'^) be the space of distributions / G iS'(R^) 
such that (1 - E Lp, with the norm 

||/||ij = ||(l-A)'=''V||„. 

Theorem 3.1. Let 1 < p < oo and k = kp = 2d\l/2 — l/p|. Then we have 
(25) MP ^ LP, Ll^ MP, 1 < p < 2 

as well as 

LP -A MP, MP ^ LP_^ 2<p<oo. 
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3.2. Canonical transformations and the associated operators. We now in¬ 
troduce a special class of operators, defined in terms of their sparsity with respect 
to Gabor wave packets and a given canonical transformation. 

A canonical transformation = xiUiV) be called tame if it enjoys the 

following two properties: 

Al. X'- is smooth, invertible, and preserves the symplectic form in R^'^, 

i.e., dx A = dy A drj] x is a symplectomorphism. 

A2. We have 

(26) \dyd^x{y,v)\<Ca,y, \a\ + \(3\ > 1, y,r] 

For example, the canonical transformation x(t, s) in 0121) satisfies Al and A2, 
with constants = Ca^y(T), provided |t — s| < T (T > 0 being arbitrary), in 
view of (IT^ . 

Let us observe that Al and A2 imply that y and x~^ globally Lipschitz. 
With the notation (z) = (H- (Japanese bracket), this property implies that 

C~\w - x{z)) < {x~^{w) -z) < C{w - x{z)) w,ze R^'^, 

for some constant C > 0 depending only on an upper bound for the first derivatives 
of X- The following class of operators was introduced in [50] . 

Definition 3.2. Let x a tame canonical transformation. Let g G iS(R‘^) \ {0}. 
We denote by FIO{x) the class of operators T : iS(R'^) —)■ such that, for 

ever^ m > 0, 

(27) \{T7i{z)g,7i{w)g)\<Cm{w-x{z))~"^, 2 :,w;eR^'^. 

We endow this space with the seminorms 

l|T||m,x= sup {w-x{z))'^\{T7r{z)g,7r{w)g)\. 

It is proved in [8] that the definition of FIO{x) does not depend on the window 
g e S{R^) \ {0}. 

The following two theorems can also be found in [8], but we provide here the 
proof for the benefit of the reader, because they were proved there in a discrete 
framework; moreover here we need some further information on the uniformity of 
the constants. 


^ In this definition we could of course take m integer; we take m real in view of possible 
extensions to the case of rough potentials, where the propagators should satisfy the same estimates 
as in (EZl) but for a fixed m > 2d. In that case it is important to allow m to be as small as possible, 
in order to reach minimal regularity. 
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Theorem 3.3. Let T G FIO{x)- Then T extends to a bounded operator on 
1 < p < oo (and in particular on Moreover, for 

every m > 2d there exists a constant C > 0 depending only on m and the dimen¬ 
sion d such that 

liril MP^MP ^ C\\T\U^. 

Proof. By the inversion formula for the STFT in fIMl) we have, if \\g\\L^ = 

T = v;v,Tv;v,. 

Since Vg : —>■ and Vf : —)■ Mp(W^), it suffices to prove 

that the operator VgTVf is bounded on Now, this in an integral operator 

in with integral kernel 

K{w,z) = {TTT{z)g,TT{w)g). 

Since T G FIO{x) we have, for every m G M, 

\VgTV;F{w)\ < \\T\U^ [ {w- x{z))-^\F{z)\dz. 

jR2d 

If we take m > 2d, the desired continuity on follows at once from Schur’s 

test and the fact that the Jacobian determinant of x is = 1, X being symplectic. □ 

Theorem 3.4. If G FIO{xi), i = 1,2, then the composition 
F/0(xi 0 X 2 )- Moreover for every m> 2d there exists a constant C > 0 depending 
only on m, the dimension d, and upper bounds for the first derivatives of Xi such 
that 

Proof. We have to estimate the integral kernel of . We write 

VgT^LTi^lV* = VgT^^\*VgT^^\*, 

so that for m > 2d 

< C||T'')||.,,„||r®||„,„ [ {w- X',(C))-’"(Xi(C) - Xi ° X 2 {z))-"'dC. 

The change of variable ( = Xi(C) ~ Xi ° X 2 {z) and the convolution property 

[ {w-C)-^{0-^dC<Cm4H-^ 

jR2d 

(valid for m > 2d, see e.g. [251 Formula (11.5)]) then give the desired estimate. □ 

The reader may want to compare the simplicity of this result with the similar com¬ 
position formula obtained in [211 Theorem A.2] for oscillatory integral operators. 
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whose proof represented really a tour de force and moreover was limited to short 
time. 

We now show that a familiar class of Fourier integral operators is in fact contained 
in FlOix) ^ suitable y. 

A real phase function <F on will be called tame if the following three properties 
are satished: 


Bl. $ e C°°(R2'^); 

B2. We have 

(28) |a“aJ<F(a;,r/)| < Q/3, |a| + |/?| >2, x, r/ G R'^; 

B3. There exists h > 0 such that 


(29) \det > S, x,r]eMf^. 

If we set 


(30) 


y = V^^{x,r]) 

i = V^<F(a:,p), 


we can solve with respect to (t, by the global inverse function theorem and 
obtain a mapping y dehned by {x,^) = x(i/, rj) and satisfying Al, A2 above as well 
as the following property: 


A3 There exists 5 > 0 such that, 

d'x ~ 

(31) |det — (i/,r/)| > 5, y,//G R'^. 

Conversely, to every transformation y satisfying Al, A2 and A3 corresponds a tame 
phase <F, uniquely determined up to a constant (see e.g. m- 
The following result was proved in [TOl Theorem 3.3]. 


Theorem 3.5. Let <F(a;, p) he a tame phase, and let x he the corresponding canon¬ 
ical transformation. Let a{x, rf) he a function in Sqq. The Fourier integral operator 

Tf{x) = (271)"'^ [ e^‘^^^’^^a{x,r])f{ri) dr] 


belongs to FIO{x)- Moreover for every m G N there exists m' G M such that 

\\T\\m,x < C||a||m' 

for some constant C depending only on m, the dimension d, upper hounds for a 
certain number of the derivatives of <F in (|28il and the lower bound constant S in 
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3.3. Semiclassical phase space methods. Here we rephrase in the semiclassical 
regime the concepts appeared above in this section. Basically we have to rescale 
function spaces and operators with the dilation operators (unitary on 

We already dehned in Introduction the spaces which are a rescaled version of 
the Sobolev spaces L^. For 1 < p < oo, consider also the modulation spaces 

Ml = {/ e : W/Wmi = \\D,-,; 2 f\\M. < +oo}. 

We have the following embeddings. 

Theorem 3.6. Let 1 < p < oo and k = kp = 2d\l/2 — l/p\. There exists a constant 
C > 0 such that 

(32) ll/lli. < ll/lUj < 1<P<2 

as well as 

(33) II/IIm-< 2<p<a). 

Proof. By fl2^ we have, for 1 < p < 2 and k = 2d\l/p — 1/2|, 

ll/llLP<qi/llM., ll/llMP<qi/llL- 

Now we replace / by Df^-i/ 2 f and we obtain ([32]), because 

IIO.-././IU. = and ||C,-.«/|Lj = 

Similarly one deduces fl3^ . □ 

Definition 3.7. Let x be a tame canonical transformation. We denote by FIOfi{x) 
the space of linear operators T : —)■ S'{W^) such that Df^-i/2TDfii/2 G 

FIO{x)? endowed with the seminorms 

ll^llm,x = \\^h-^/^TDf^i/ 2 \\m,x- 

From Theorems 13.31 we obtain at once the following result. 

Theorem 3.8. Let T G FIOfi{x). Then T extends to a bounded operator on 
1 < p < oo. Moreover, for every m > 2d there exists a constant C > 0 
depending only on m and the dimension d such that 

WTWmi^mi < C\\T\\l^^. 

Corollary 3.9. Let T G FlOnix), I < p < oo and k = kp = 2d\l/p — l/2\. Then 
T extends to a bounded operator T : —)■ z/ 1 < p < 2 and T -. L^ ^ for 

2 < p < oo. Moreover for every m > 2d there exists a constant C > 0 depending 
on m, d, p such that 


||T||iP^^,<q|T||^,^ forl<p<2 
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and 

< CIITII^,^ for2<p<oo. 

Proof. The results follow at once from Theorem 13. 8l and Theorem 13.61 For example, 
for 1 < p < 2 we have 

< c,C2C3||r||»_j|/||i.. 

□ 


The following result will be crucial in the sequel. 

Theorem 3.10. //Tb) e FlOniXi), '^ = I 52 , then the composition is in 

FIOti{xi°X 2 )- Moreover for every m > 2d there exists a constant C > 0 depending 
only on m, the dimension d, and upper bounds for the first derivatives of Xi such 
that 

Proof. We write 

= Df,-^/2T^^^Dni/2Dn-i/2T^^^Dni/2 

and apply Theorem 13.41 □ 

4. The Schrodinger flow: wave packet analysis 

In this section we prove that the oscillatory integral operators appearing in Sec¬ 
tion 2 belong to the class FlOnix^it^ ■s)) for a suitable canonical transformation x^ 
dehned below. Then we will use Theorem 13.101 to treat their composition for large 
time. 

We adopt the notation of Section 2. In particular the oscillatory integral opera¬ 
tors are well dehned for |t — s| < 6. 

Proposition 4.1. Let |t — s| <6 and a G Sqq. Consider the operator 

Then T can be written as a Fourier integral operators 

(34) Tf{x) = (271)-" / f, 5 , a:, hp)/(p) dp 

Jm.'i 


where 


^{t,s,x,p) = yp + S{t,s,x,y) 
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is the generating function in the coordinate^ x, rj, and a E Sqq. Moreover for every 
m G N there exists m' eN such that 

||6(/i, S, •, ■)||m ^ C*!!Cl||m' 

for some constant C > 0, depending only on m, the dimension d, upper bounds 
for a certain number of the derivatives of (t — s)S(t, s, x, y) in (fT^ and the lower 
bound constant S in m- 

Proof. Set S(t,s,x,y) = (t — s)S(t, s,x,y) and fj = h(t — 3 ) 7 ]. By the Plancherel 
theorem we can write T in the form 


Tf(x) 

- . .. ^ f { [ exp(i(ft(i-s)r‘|OT + S((,s,x,j,)])a(x,j,)d»)7(„)rf-). 

{2Tii{t-s)hp^ \J^d V / / 


We now apply the stationary phase principle in the form proved in page 320 
and Lemma 3.2]; the assnmptions are satished becanse for |f — s| < 5 we have, by 
(m and ([ED, 

\dfd^S{t,s,x,y)\ <Ca,p |a| + |/3| > 2, 


and 


det 


d^S{t,s,x,y) 

dy'^ 


> 5>t). 


Moreover, since y = y{t, s, x, rj) is the nniqne solntion of —rj = dS{t, s, x, y)/dy, the 
fnnction y = y{t, s, x, f}/(t—s)) will be the nniqne solntion to —fj = dS(t, s, x, y)/dy, 
and we obtain 


= exp (^i{h{t - s))~^[y{t, s, x, fj/{t - s))fj + S{t, s, x, y{t, s, x, r//(t - s)))] j 

X b{h, t, s, X, f])) 

for some amplitnde b{h,t, s, ■, ■) belonging to some bonnded snbset of Sqq, when 
\t — s\ < S, 0 < h < 1. Again, we can write the last expression as 

exp (^ih~^[y{t, s, x, hri)hri + S{t, s, x, y{t, s, x, hp))]^b(h, t, s, x, (f — s)hr]), 

which gives the desired expression for Tf. 

The estimates of the seminorms of b in terms of those of a follow from the proof 
of the stationary phase principle. □ 


^Namely, here y = y(t, s, x, g) is the unique solution to —y = dS(t, s, x, y)/dy. 
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Corollary 4.2. With the notation of Proposition 
with 


we have T G FIOn{x^{t, s)), 


s){y, v) = h s, hfl‘^y). 


Moreover, for every m G N there exists m' G N such that, for — s| < 5, 

nit,< cii«iU' 

for some constant C > 0 independent of a and h, t,s (0 <h<l,t — s<6). 

Proof. We have to prove that Df^-i/2TDf,i/2 G FIO{x^{t,s)). Using (IM|) and a 
scaling argument we can write 

D^-,i2TD^,i2f{x) = ( 27 r)-'^ [ t, s, h^^^y)f{y) dy. 


Now, the phase h~^^(t, s, generates the canonical transformation x^(t, s) 

as in the statement. Moreover h~^^{t, s,h}/‘^x,h}/‘^y) is tame uniformly with re¬ 
spect to h and t, s, for 0 < h < 1 and |t — s| < S, in the sense that the re¬ 
quired bounds hold with constants independent of these parameters, or equiva¬ 
lently x^(t, s) satishes the properties Al,A2,A3 in Section 3.2 above, with uniform 
bounds. This is clear, because x(t, s) satishes Al,A2,A3 uniformly with respect 
to t, s for |t — s| < 6 (possibly for a smaller value of 5): A1,A2 follow from ([13]), 
whereas A3 holds (even in the presence of a magnetic held as in Section 6 below) 
by |52l Proposition 2.3’] applied with a = (3 = t). 

Finally, the symbol b{h,t, s,h}^'^x,h}^‘^y) has seminorms in S'qq dominated by 
those of b{h, t, s, •, •). Hence the desired result follows from Theorem 13.51 □ 


Remark 4.3. With the notation of the previous corollary, notice that for any 
s,T,t eM. we have 

X^{'t,s) = x^(f,r) o/(r,s). 

Indeed, this follows at once from the case h = 1 and the definition of x^- 

5. Proof of the main result (Theorem 11.11) 

Let us hrst prove ((8|) and ((9|). We know from Section 2 that the propagator U{t, s) 
is an oscillatory integral operator of the form fl^TD for |t — s| < 6 , and therefore by 
Corollary 14.21 we have U{t,s) G FIOfi{x^{t,s)), with seminorms \\U{t,s)\\^^n(^i-gy 
m G N, uniformly bounded with respect to h, t, s for \t — s\ < 6. Using the evolution 
properties of the propagator. Remark 14.31 and Theorem 13.101 we see that, for any 
T > 0 we have U{t,s) G FlOnix^it, ^)) with seminorms 

(35) 

for a constant Cq independent of h, t, s, for \t — s\ <T (but depending on T); here 
we are using the fact that constant C in Theorem 13.101 can be chosen independent 
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of such parameters, being -s) uniformly tame. We then deduce (0) and (l9i) 
from Corollary I3.91 

We now prove (1T0|) and (ITT]) . By arguing as above it suffices to prove that 

(36) s C'(r)«"u)(n)''+‘|( - s| 
for |f — s| < T and for some m > 2d. 

Now, by Corollary 14.21 and flT^ the operator G^^\t,s) defined in fflSj) belongs 
to FIOfi{x^{t, s)), with seminorms 

(37) - si"*' 

for a constant Cm independent of h, t, s, for \t — s\ < 6. On the other hand we have 

(t, s)f:= (t, s)f- U{t, s)f = -ih-^ U{t, (r, s)fdT 

so that by (l35|) . (1^ and Theorem 13.101 we obtain 

(38) s) t,,,.;,,,) < 

We then can write 

(O, t, s)-U{t, s) = {U{t, tL-i)) ... {U{C, s)+R^^\c, s))-U{t, s) 

and argue as in Lemma 3.2], applying the composition estimate in Theorem 
13.101 We sketch the argument for the benefit of the reader. 

One expands the above product and obtains a sum of ordered products of opera¬ 
tors, where each product has the following structure: from right to left we have, say, 
qi factors of type U, pi factors of type R^^\ q2 factors of type U, p2 factors of type 
R^^\ etc., up to qk factors of type U, pk factors of type R^^\ to finish with qk+i 
factors of type U. Here pi, ■ ■ ■ ,Pk, di, ■ ■ ■ dk, qk+i are non negative integers whose 
sum is L, with pj > 0 and we can of course group together the consecutive factors 
of type U, using the evolution property of the propagator. Now, we estimate the 
seminorm || • \\mx'^{ts) ordered product, for some fixed m > 2d, using 

Theorem 13.101 combined with the known estimates for each factor, namely fl35l) and 
fl38|l . which we rephrase here as 

(39) < Coh^\tj - 

If C denotes the constant in Theorem 13.101 we then dominate the seminorm || ■ 
Wmx’^its) above ordered product by 

k Pj 

coh^\tJ^+^ - 

j=i i=i 
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where Jj = pi + .. .+pj_i + qi +.. . + qj for j >2 and Ji = qi. Now, it was proved in 
the last part of the proof of [211 Lemma 3.2] that the snm over pi,... ,pk,qi,..., 
of these terms is in tnrn dominated by 

c{T)h^uj{n)^+^\t - s\ 

for |f — s| < T. This gives fl36|) and conclndes the proof. 


6. The case of magnetic fields 

In this section we extend some of the previons resnlts in the presence of a mag¬ 
netic held. Hence, consider the Schrodinger eqnation 

1 2 

(40) ihdtu = 2 ( ~ ~ u + V{t, x)u 

where V{t, x) and A{t^ x) = (Hi(t, x),..., Ad{t, x)) are electric scalar and magnetic 
vector potential of the held, t G M, x G 

Assnme the following hypothesis (cf. [2411^ 152]). 

Assumption (B) 

a) For j = 1,... ,d, Aj(t, x) is a real function of (t, x) G M x and dfAj(t, x) 
is in {t, x) G M X for every a G N'^. Moreover there exists e > 0 such 
that 

\d^B{t,x)\ < Ca{l + \x\)~^~'', |a| > 1 

\dfA{t,x) \ + \d^dtA{t,x)\ < |a| > 1, 

for (t,x) G M X where B(t,x) is the magnetic field, i.e. the skew- 
symmetric matrix with entries Bj^k(t,x) = {dAk/dxj — dAj/dxk)(t,x). 

b) V{t,x) is a real function of (t, x) G M x with dfV{t,x) continuous in 
{t, x) G M X W^, for every a G and satisfying 

\dfV{t,x)\ < Co,, |a| > 2, (t,x) G M X R'^. 

It was proved in |52l Sections 2,3] that if H = 0, nnder Assnmption (B), all the 
short time resnlts abont the propagator, as well as the constrnction of the operators 
E^^\t,s), G^^\t,s) snmmarized in Section 2 continne to hold if A > 1. To be 
precise, the same formnlas for E^^\t,s), G^^'^{t,s) and E^^\kl,t, s) hold, where 
now the amplitndes Ok satisfy the transport eqnations 


(41) 


dak 

dt 


^ 


X 


J=1 


dak 

dx.i 


+ ^(^A,,S{t,s,x,y) - - div^j, A(t,x)jafc = ^A^-afc-i, 
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with ao(t, s, x,y) = 0 and the initial conditions ai(s, s, x,y) = 1 and ak{s, s, x, y) = 
0 for k = 1 , 2 ,... (these transport eqnations hold even when the electric potential 
is present!!). The key estimates (IT7|) - (l2^ are valid too, for iV > 1. 

We have therefore the following resnlt. 

Theorem 6.1. Under the above Assumption (B), the conclusions of Theorem \l.l\ 
still hold, at least for N > 1. 


Proof. As we have already observed, the resnlts in Section 2 keep valid nnder the 
Assnmption (B), if 1/ = 0 and A^ > 1, therefore the whole snbsequent analysis 
applies, giving the desired result if Id = 0. 

Consider now the case when a potential V{t,x) is present. Let us prove just 
part b) of Theorem 11.11 (the hrst part is similar and easier). We apply the gauge 
transformation 

Q{t)u{t,x) = eyi\){ih~^G{t,x))u{t,x) 

which maps the solution of the equation with potentials A, V into the solution of 
the equation with potentials A + V — dtG, cf. [S2]. More precisely, if we 

denote by U(t, s) and U'(t, s) the propagators, with the prime denoting that of the 
new equation, we have 

u{t,s) = g{t)-^u'it,s)g{s). 

In particular if we choose G{t, x) = V(T , x) dr, in the new equation the electric 
potential is absent and and magnetic potential A' = A + dxV{T, x)dT, still 
satisfying the above Assumption (B). Since the Lagrangian function L{x,v,t) = 
+ A(t, x) ■ V — V(t, x) changes to L' = L + dG(t, x(t))/dt, the flow remains the 
same and the new generating function is 

V (r, x) dr — f V (r, y) dr 

Jo 

where S{t,s,x,y) is the generating function for the original equation. The am¬ 
plitudes Ofc for the new equation are the same as those for the original one, be¬ 
cause the transport equations ga do not change under the substitution A —)■ 
A + Jg Va;ld(r, x) dr, S ^ S', V ^ 0. As a consequence we have 

E^^\t,s) = g-\t)E^^y{t,s)g{s) 


S’{t,s,x,y) = S{t,s,x,y) + 


^Indeed, since S satisfies the Hamilton-Jacobi equation dtS + i(Va;5' — A)^ + V = 0, we have 


fhdt - -{-ihVx - -V{t,x 


{2'Kih{t — s))^/^ 




(2TTih{t — s))‘^P 


Jh ^S(t,s,x,y)^j^ 


B'i 


dta+{VxS-A)-Vj;a+^{A,^S--J^ -div^, d)-^A 3 ;aj dy. 

2 t — s 2 J 
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and therefore 

= g-\t)E^^y{n,t,s)g{s) 

where the the primes denote the corresponding operators for the new eqnation. As 
a conseqnence 

E^^\n, t, s) - u{t, s) = g-\t) s) - u\t, s))g(s). 

Now, we know that fl36l) holds for the difference E^^^' s) — U'(t,s) and we 
want to prove that the same is trne for E^^\il,t, s) — U{t,s). Using Theorem 
13.101 it is then snfficient to prove that g{t) G EIO^^x^) for some tame canonical 
transformation x- This follows by Corollary 14.21 becanse g{t) can be written in 
the form fIMD with <h(t, x, 7]) = xrj + V(r, x) dr, 6=1, and 

X(v,r,)={v,r, + j\,V(T,x) dr). 

□ 

Remark 6.2. lUe do not know whether the conclusion of Theorem \1.1\ holds for 
= 0 m the presence of a magnetic field as above (in fact, in /5li Theorem 5] 
the case N = 0 is excluded as well). Indeed, the approach in 1^7] / requires very 
precise estimates on the generating function S{t, s,x,y), which seem to he quite 
tricky when a magnetic field is present. 


7. Concluding remarks 


7.1. Sharpness of the results. In general we cannot avoid the dichotomy p > 2, 
p < 2 in Theorem ll.il Consider, for example, the case of the harmonic oscillator. 


(42) 


1 « 1 , ,2 

lOiU = —Am H— \x\ u 
2 2 ' ' 


(we take h = 1 for simplicity, so that = Lf(f). By the Mehler formnla (see 
e.g. |3Sl formnla (5a)]), the propagator U{t, 0) at t = 7r/2 is (np to constant factor) 
jnst the Fonrier transform / eA /. Now, we have the following resnlt. 


Proposition 7.1. Let 1 < p < oo, ki,k 2 G M. The Fourier transform maps 
—)■ continuously if and only if 

ki > 2d{l/p — 1/2) and /c 2 < 0, if 1 < p < 2 


/ci > 0 and /c 2 < —2d{l/2 — 1/p), if 2 < p < cx). 
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Proof. Sufficient conditions. The desired estimates follow by the inclusion relations 
of Sobolev spaces and duality if we prove the continuity of the Fourier transform 
—)■ U' with k = 2d{l/p — 1/2), for 1 < p < 2. This is equivalent to the estimate 

ii(i+i?r)-'‘'TK)iiL.<cii/iiL.. 

which in turn is a consequence of the Hardy-Littlewood-Paley inequality [3l The¬ 
orem 1.4.1] (that is the same estimate with (1 -|- replaced by in the 

left-hand side). 

Of course the result also follows from Theorem 11.11 applied to the equation fl42|) . 
with s = 0, t = 7r/2 {h = 1). 

Necessary conditions. Let us first prove the condition on k 2 when 2 < p < oo. 
By duality this will give the condition on ki for 1 < p < 2 as well. 

Consider the space 

A={f estm"): /K) = 0for|J|>l}. 

By Bernstein inequalities (see e.g. |53l Proposition 5.3]) there exists a constant 
C > 0 such that 

II/IIl- <c\\fh. v/gA 

fci 

even for fci < 0 (for fci > 0 this is trivially true for every / G iS(M'^)). 

Hence, if the Fourier transform is bounded —)■ we have 

(43) <C||/|Up V/gA 

Now, suppose by contradiction that ^2 > —2d{l/2 — 1/p). Fix / G {0} and test 

this estimate on f{x/X), with A > 1, so that f{x/X) belongs to A too. We have 

(44) IK1-A)‘>/=7(7a)||„<C||/(./A)U. 
where 

(1 _ = A'^+^2[(A-2 - A)’^^/^f]{X-). 

By dominated convergence theorem {\x\^^ is locally integrable, because ^2 > —d) 
we have [(A“^ — f]{x) —?■ \D\'^^f{x) for every x G Mf, and therefore by the 

Fatou theorem, 

0 # ll|C|‘“/l|i. < liminf ||(A-= - A)‘-/"/||„. 

A^+oo 

Hence, letting A —>■ -l-cxo in (1441) we obtain k 2 < —2d{l/2 — 1/p), which is a contra¬ 
diction. 

Let us now prove that /c2<0 ifl<p<2 . This will give the condition on ki for 
2 < p < cx) as well. 
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Fix / G \ {0}, and let fx{x) = f{xi — X,X 2 , ■ ■ ■ ,Xd), A > 0. Suppose the 

following estimate holds: 

(45) \\MLl^<C\\fx\\Li^=C\\fh.. 

We have 

(1 - = (1 - A)"^/2[e-*"«i/(0] 

= [ e“^(l + + A),-X2,...,-Xd) da; 

jRd 

= [ e-*"«(A-2 + |1 + X-^yi\^ + + • • • + f{y) dy. 

As A —>■ +CX) the last integral converges to f{^), for every ^ G by the dominated 
convergence theorem. One can then conclude as above, by letting A —)■ +oo in 0451) 
and using Fatou theorem. 

□ 

7.2. The case of rough potentials. Here we present an extension of the hrst 
part of Theorem II.II to a class of potentials in the so-called Sjostrand class 

^oo,i(^d) ^ ||/|Uoo,i(M.) := [ \\VJi;0\\L^d^ < +oo} 

jRd- 

where, as usual, g G \ {0}. 

Functions in this space are bounded in and locally have the mild regularity of 
a function whose Fourier transform is in L^. Hence, it is not required the existence 
of any derivative. We refer to the books [T51I25] for a detailed analysis of this 
function space and its role in Time-frequency Analysis and Mathematical Physics. 
Here we just need the notion of narrow continuity for a one-parameter family t ^ ft 
of functions in this space. In fact, it is quite a weak notion of continuity. 

Definition 7.2. Let / C M 6e an interval. We say that a map I 3 t ft E 

is continuous for the narrow convergence if it is continuous in 
(weakly) and if there exists a function H G such that sup^jg^d \Vgft{^^0\< 

dd{0 for every t ^ I and almost every ^ G 

Consider now the equation (with h = 1, for simplicity) 

idtu = (—A)''^^n -|- I4(t, x)u 3- Viit, x)u + Vo{t, x)u 

with 0 < n < 2, t G M, a; G Suppose 

•) e for |«| = j = 0,1, 2, 
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with V 2 and Vi real-valued. Assume moreover that the map R 9 t ha d^Vj{t, ■) G 
M°°’^(R'^) for |a| = j, j = 0, 1 , 2, are continuous in the sense of the narrow conver¬ 
gence. 

Under this assumption the propagator U{t,s) was constructed in [13] (see also 
[ 9 IITT] and |32| for a similar class of potentials having Fourier transform in L^) 
and shown to be bounded — )■ for every 1 < p < cxo. As a consequence of 
Theorem 13.11 we then see that U{t, s) enjoys the continuity property 

U{t, s):Ll^LP, 1 < p < 2, U{t, s) : ^ 2 < p < 00 

with k = kp = 2d\l/2 — l/p|. 
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